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Abstract 

In this paper we compute the three-loop corrections to the (5 function in a mo- 
mentum subtraction (MOM) scheme with a massive quark. The calculation is 
performed in the background field formalism applying asymptotic expansions for 
small and large momenta. Special emphasis is devoted to the relation between the 
coupling constant in the MOM and MS schemes as well as their ability to describe 
the phenomenon of decoupling. 

It is demonstrated by an explicit comparison that the MS scheme can be con- 
sistently used to relate the values of the MOM-scheme strong-coupling constant 
in the energy regions higher and lower than the massive-quark production thresh- 
old. This procedure obviates the necessity to know the full mass dependence of the 
MOM (3 function and clearly demonstrates the equivalence of both schemes for the 
description of physics outside the threshold region. 



PACS numbers: 12.38.-t, 12.38.Bx, 14.65.-q 



1 Introduction 



Within the perturbative framework, the MS scheme PQ [2] based on dimensional regulariza- 
tion [31 m [3] is a well-estabhshed scheme for the renormahzation of fields and parameters. 
This applies in particular to a^, the coupling constant of Quantum Chromodynamics 
(QCD). One of the major advantages of the MS scheme is its simplicity in practical 
applications. The main reason for this is that it belongs to the class of so-called mass- 
independent schemes where the renormahzation constants are independent of the precise 
configuration of masses and external momenta involved in the problem. 

Within the MS scheme, the beta function governing the running of as is know in 
the four- loop approximation [6], [7]. In order to correctly account for the heavy-quark 
thresholds, also the corresponding matching (or decoupling) conditions are needed, which 
allows for a precise relation of at widely separated energy scales like, e.g., the tau lepton 
and Z boson masses. Four-loop running goes along with three-loop matching, which is 
also known since more than ten years [S]o 

Other renormahzation schemes which do not have the nice property of mass- 
independence are significantly more complicated from the technical point of view — 
mainly because one has to deal with Feynman integrals involving many mass scales. Still, 
at the level of precision which has been reached in the recent years, it is necessary to have 
a cross check of the dependence on the renormahzation scheme. In this paper, we want 
to provide an alternative set-up to the running and decoupling of as in the MS scheme 
and consider a momentum subtraction (MOM) scheme for the definition of as- We will 
provide MS to MOM conversion formulae and the MOM beta function in the three-loop 
order and are thus able to cross check the MS running of A two-loop analysis has been 
performed in Ref . [H] . In this paper we check the calculation of Ref . [H] and extend the 
analysis to three loops. 

The remainder of the paper is organized as follows: In the next section, we describe 
our setup. In particular, we derive the relation between the strong coupling in the MS 
scheme and in two versions of the momentum subtraction scheme and provide the cor- 
responding beta functions. In Section [21 we present our analytical results for the gluon 
polarization function in the background field formalism and discuss the phenomenological 
applications in Section HI where we compare the running in the MS and MOM schemes. 
Our conclusions are summarized in Section [51 

2 The strong coupling in the MOM scheme 

For convenience, we adopt Landau gauge, which has the advantage that the renormahza- 
tion group equations for the gauge parameter and as decouple. Furthermore, we require 
that the polarization function of the gluon vanishes for = —q^ = > 0. 

For the practical calculation, we adopt the background field gauge [12], which has 
the nice feature that the (3 function of the strong coupling is determined from the gluon 

^Recently, also the four-loop decoupling constants have been computed [51 llOj. 
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polarization function alone. The latter is given by 

n^^(g) = (-^?^V + gV) n(g2) , (i) 

which is conveniently decomposed as follows 



i>l 



In the z-loop contribution, the dependence on q, fi and the various quark masses is ex- 
plicitly displayed. Formulae ([1]) and ([2]) hold both in the MS and MOM schemes. The 
corresponding functions, n(g^) and n'^*-'^(g^), can be used to obtain a relation between 
as and a^^^, the strong couplings in the MS and MOM schemes, using the fundamental 
concept of the invariant charge 



1 + n^ioM(^2) I + n(g2) • 

It is an important and unique feature of the background field gauge that the invariant 
charge is expressible in terms of the coupling constant and the gluon polarization op- 
erator only in exactly the same simple way as in QED. We define U^^^[q'^) such that 
jjMOM^_^2^ _ g g^Yid, Consequently, we have 



1 + Ci hC2 +C3 



C2 = -n;,^^+ n 



(2) , /'qW'^ 



C3 = -n?) + 2n«nf)- n(^)) , (4) 
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where Ilg^ = 11*^*) (—/x^) has been introduced. It is instructive to look at the explicit 
expressions in the massless limit with Uf = ui massless quarks. In this case, we obtain 



«^V^2W205 10 \ /a^)(/i2)\ /90391 
1 : -tt; TTiT-i 
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20423^,^, 2320^, / 
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(5) 
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In analogy to the MS scheme, the f3 function in the MOM scheme is defined through 



d a 



MOM 



d/x2 



71 



MOM/ MOM\ 



a 



MOM\ 2 



TT 



MOM / "^s 



i>0 



MOM 



TT 



(6) 



where — in contrast to the MS scheme — the coefficients Pf^'^^ are functions of the 
renormahzation scale /i and the quark masses. With the help of Eq. ([3]), where we replace 
on the right-hand side the MS renormalized quantities by the bare ones, it is possible to 
obtain a relation between and the coefficients n^^)'^^^ = U^''^'^^^{q\ jj,^, {M^}), 

which reads 



/3 



MOMj^^MOM-j 



a 



MOM 



MOM 



TT 



i-E.,.(.-i)(^)'n 



(i),MOM 



(7) 



From this equation, one can easily derive convenient formulae for f3f^'^^. Note that the 
term in the denominator of Eq. (jTj) contributes for the first time at the three-loop order. 
Let us also mention that, starting at this order, a non-trivial dependence occurs on the 
right-hand side of Eq. ([7]) which has to cancel in the proper combination of the ]jMMOM 
functions. 

The functions and are known analytically [TT]. The three-loop contri- 

bution is evaluated in the asymptotic regions for large and small quark masses 

analytically in this paper. An approximate formula valid for arbitrary quark masses is 
easily obtained by interpolation between the low- and high-energy regions. 

In the massless limit, the first three coefficients are given by 



nMOM 
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6 

32 



C(3) C\Tni 



+ ^C(3) CpT'Tif + 2Cj,Tni 



C(3) ) CaT'u] 



where Ca = 3,Cir = 4/3, T = 1/2 and rii is the number of massless quarks. Since 
the first two coefficients of the /? function are scheme independent P^^i^ and 



coincide with their counterparts in the MS scheme. 0^^\^-, however, differs from its MS 
counterpart [TJl |TB]. It is worthwhile to mention that 0^^\^ contains the Riemann C, 
function C(3), which in the MS scheme only appears at the four- loop order. 
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The three-loop results in Eqs. ([5]) and ([8]) are new, and the two-loop expressions are 
in agreement with Ref. [llj . 

The practical evaluation of 11^°^ (g^) entering the equation for the beta function can 
be reduced to the evaluation of n(g^) in the MS scheme. The corresponding relation is 
obtained from Eq. ([3]), this time for arbitrary values of and /x^, which can be solved for 
■qmom_ p^^iqy properly replacing q;^°^ by using Eq. (jl]), one gets (the dependence on 
the quark masses is suppressed) 



]-[(l),MOM(^^2^| 
J-[(2),MOM^^^2^| 
]-[(3),MOM(^^2^) 



n(2)(g2 

n(3)(g2 



n, 



-n, 



(1) 



(2) 



(9) 



where Ilg^ is defined below Eq. (jlj). Note, that by construction we have n'^"^(— /x^) = 0. 

The polarization function in the MS scheme is obtained in the standard way by renor- 
malizing in the MS scheme, the quark masses in the on-shell scheme and taking care 
of the gluon wave function renormalization. 

In Ref. p^, it has been observed that there are relatively large coefficients in the 
relation between and when running from as{Mz) down to, say, ^^(Mt-). The 

situation was improved in Ref. [llj by a simple trick of rescaling the scale parameter /i. 
Let us start from the massless limit corresponding to fi ^ Mf. In this case, relation 
assumes the form 



1 + 10.417 ' ^ + 126.350 ' 



An 



(6)^,2^ 



+ 2000.062 



as 



3 n 



(10) 



In a next step, we introduce a new, rescaled MOM scheme with the help of 



or, equivalently (with L = In(xQ)), 
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where Pf^^i^ are given in Eq. ([8]). The corresponding generahzation of Eq. (jlj) reads: 



1 + Ci hC2 +C3' 



-02 = -n(n(n«) -2nWr, + r„ 

C3 = -n(n2n«nf-(n«)V3(n«)%, 

-2nfVi-3ni'^2 + r3. (13) 

In a next step, following Ref. [11], we tune the parameter xq so that the difference between 
^MOM ^(6) gi^g^j^^g only in order a^. The result read^ 

125 

\n{xl) = — , xo^ 2.1044, (14) 

which leads to 

ri ^ -2.60417, 
r2 ^ 4.3635, 

rs ^ 2.2313. (15) 

It is instructive to look again at the relation between a^*-"^ and ai^^ for /i ^ Mj, 
which is now given by 



2 / „ \ 3 



where 

t = ™55_5ZZc(3) « -20.8472, 

168 8 ^ ^ 

~ 101389 345779,, , 222305,, , , , 

As compared to Eq. ( fTOl) . we observe a significant reduction in the magnitude of the coeffi- 
cients in the rescaled relation m, both at the two- and three-loop ordersl Furthermore, 



^Note that there seems to be a misprint in the nmnerical value of xq quoted in Ref. [llj . however, in 
the caption of Fig. 4 therein it is correct. 

•^Note that our value for the two-loop coefficient in Eq. (fT6| (—20.8472) differs from the one obtained 
in Ref. [II] (-32.46). 
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there is a different sign in the three-loop coefficient as compared to the two-loop one, 
which points to a better convergence of the perturbative expansion. 

In the massless limit, both definitions for Q;^°^(/i^), namely Eqs. ( fTTl) and (fT2l) . are 
completely equivalent. Following again Ref. [H], we choose Eq. ( IT2l) as the proper defini- 
tion of the for a// values of /i. This choice has the advantage that the thresholds 
in the corresponding function /J'^om j-emain "physical", that is located at — /z^ = 4Mq. 
This follows directly from the relation between both j3 functions: 



nMOM _ oMOM 
Po ~ Po ' 

^MOM _ ^MOM 



It is interesting to remark that the rescaling procedure significantly improves the MOM 
to MS relations also for moderate and even rather low values of (see below). 

In the applications of Section HI we consider the strong coupling both for energy scales 
of the order of or larger than and for those significantly smaller than the top-quark mass. 
In the latter case, we construct different MOM and MOM schemes, which are derived 
from the choice rij = 5 as the massless limit. In this case, we obtain the values 



Xo = 


g415/552 _ 2.1208, 


k = 


140689 1699,,,, 

C(3) ^ - 

1656 24 ^ 


h = 


143409281 1225793 
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415 

-2.8819, 


r2 = 


4.6719. 

476928 


r3 = 


1188703175 70508^ 
68677632 ' 55296 



-0.1385, 

518435 



576 



C(5) ^ 831.5896 



■C(3) ^ -1.9809. (19) 



The same comments and conclusions hold as for Uf = 6. 



3 Results 

Let us in a first step briefly describe the evaluation of the gluon polarization function up 
to three loops within the background field formalism involving heavy quarks with generic 
mass Mq. The basic idea is to evaluate n(g^) for large and small external momenta and 
to obtain an approximation for all values of /Mq by a simple interpolation procedure. 
Note that in our case the external momentum is space-like so that there are no problems 
with particle thresholds. Up to the two-loop order, only one quark flavour can occur in a 
diagram. At three loops, there are diagrams with a second closed fermion loop so that in 
principle a further mass scale can occur (see, e.g., the diagram in Fig. 1(g)). However, we 



6 




Figure 1: Sample diagrams contributing to the gluon propagator in the background field 
formalism at the two- and three-loop orders. Diagrams (a), (b), (e) and (f) only contain 
massless lines, while the others also contain massive ones due to the presence of the 
heavy-quark loop (thick line). (We have used the package JaxoDraw [TTl [H] to draw the 
diagrams.) 

assume a strong hierarchy in the quark masses such that we can always neglect the lighter 
mass. Thus, in this section, we consider QCD with total number rif of quark fiavours. 
One quark, Q, has the (pole) mass Mq, and all other ni = nj — 1 quarks are considered 
as massless. 

As mentioned above, the MS renormalized polarization function is needed in Landau 
gauge. However, in our calculation we adopt a general gauge parameter ^ since the 
complexity is comparable to Landau gauge. 

Some sample diagrams for n(g^) are shown in Fig. [T] The diagrams are divided into 
two classes: completely massless diagrams and and diagrams involving massive-quark 
loops. The only scale in the massless diagrams is the external momentum. Thus they can 
be evaluated using MINCER [191 [20]. Iii the second class, the mass of the heavy quark sets 
another scale which makes the calculation significantly more difficult. The one- and two- 
loop calculations can be performed analytically, and the results can be found in Ref. pT] . 
At the three-loop order, however, an exact calculation is not yet possible. We perform 
an asymptotic expansion in the limits <^ Mq and ^ Mq. Note that due to the 
diagrams containing massive-quark loops along with massless cuts (see, e.g.. Figs. [D^d) 
and (h)) also the small-g^ expansion turns out to be nontrivial. As a result, one encounters 
ln(g^/MQ) terms also in this limit. 

All Feynman diagrams are generated with QGRAF [2T] . The various diagram topologies 
are identified and transformed to FORM [22] with the help of q2e and exp [23t [2l] . The 
program exp is also used in order to apply the asymptotic expansion (see, e.g., Ref. [25] ) 
in the various mass hierarchies. The actual evaluation of the integrals is performed with 
the packages MAT AD [26] and MINCER [20], resulting in an expansion in d — 4 for each 
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diagram, where d is the space-time dimension. 

We computed four expansion terms for small and six terms for large external momen- 
tum. In the following, we present only the leading and subleading terms of the corre- 
sponding expansions for the gluon polarization operator and the MOM jS function for 
yU^ = and provide the complete expressions in a Mathematica fileEl For completeness, 
we also list the one- and two-loop results, which agree with the corresponding expansions 
of the exact expressions [11]. It is convenient to cast the result in the form 



^(g^) = ^-l(g^) + ^-^(g^M^ 



(20) 



and introduce the variable 



(21) 



The results for the massless MS renormalized polarization function reads 
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See http : //www-ttp . particle . uni-karlsruhe . de/Progdata/ttp08/ ttp08-50 
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and in the large-Z region, we obtain 
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The small- and large-Z expansions of the MOM /3 function read: 
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Note that our result for the MOM f3 function exphcitly demonstrates the vahdity of 
the Applelquist-Carazonne theorem [27] at the three-loop level. Indeed, one can easily 
check that, for i = 0,1 and 2, one has 



lim 3, 
lim (3. 

Z^oo 



MOM 



/oMOM/ 



MOM 



MOM, 
i,ml > 



ni), 
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where ('^z) the three-loop contribution to the MOM (3 function in the massless 

limit (see Eq. ®). 

In Fig. [2], we present the results for the MOM j3 function in graphical form, where 
the one-, two- and three-loop coefficients are shown as functions of Z in the Euclidian 
region. Next to the low- and high-energy approximations (dashes) including the Z^ and 
1/Z^ terms, also the interpolation functions (dotted) are shown. At the one- and two- loop 
orders, these results are compared against the exact result (solid line). For demonstration 
purpose, we have chosen n/ = 5 at the three-loop order. Very similar results are obtained 
for other values of rii. 



4 Phenomenological applications 

In the following, we discuss the numerical impact of the results obtained in this paper. 
In particular, we consider the MS quantity af\Mz) as input value and evaluate the 
strong coupling at lower and higher energy scales with different numbers of active flavours. 
On the one hand, this can be done in the MS scheme applying the usual running and 
decoupling procedure (see, e.g., Refs. [H [28]). In this case, one has to specify a scale /ig 
where the heavy quark Q is integrated out. On the other hand, it is possible to switch 
from the MS to the MOM (MOM) scheme for fi = Mz and perform the running with the 
help of the MOM (MOM) (3 function. The results obtained at lower and higher energies 
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Figure 2: The difference (3f^^^{Z) - pf^'^^iO) {i = 0,1,2) as a function of Z in the 
Euchdian region. In the three-loop resuh, ni = nf — 1 = 5 has been chosen. In each 
frame, the dashed lines represent the approximation for small and large values of Z, the 
dotted curve is the result of the interpolation, and the solid line (for and is 

the exact result from Ref. [TT] . 
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Table 1: and a^^^ for various values of fi from region A. As input, as{Mz) = 0.118 
is used, which is transformed with the help of Eq. (jl]) to a^°^(M^). The running to 
different values of /i is achieved with the help of the appropriate /3 function. a^°^(Q;s(/i)) 
is obtained from ^^(/i) using Eq. (jl]). 
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Table 2: as and a^'^^ for various values of /i from region A. As input, as{Mz) = 0.118 
is used, which is transformed with the help of Eq. ( fT3l) to q;^°^(M^). The running to 
different values of is achieved with the help of the appropriate (3 function. a^^^{as{fi)) 
is obtained from ^^(/i) using Eq. (IT^ . 



can also be translated back to the MS scheme, and a comparison can be performed. In 
this way, we can check the consistency between the two renormalization schemes. 
For our numerical analysis, we use the following input values 

a(5)(M^) = 0.118, Mb = 4.7 GeV, Mt = 175GeV, (27) 

where Mq represent the pole quark masses and for the decoupling scales we choose fig = 
2 Mq. The running and decoupling in the MS scheme is performed with the help of 
RunDec ^29j . 

We consider two regions of energies. Region A starts from the Z-boson mass, Mz, 
and extends to energies much higher than the top-quark mass, say, 1000 GeV. In this 
region, we investigate the evolution of the strong-coupling constant in the MOM and 
MOM schemes with five massless and one heavy quark, the top quark. Thus in total six 
quarks are present in the theory which we denote as n^^^^ = 6. On the other hand, in 

region B, we consider the evolution of a^^^(/i) and a^°^(/i) from /i = Mz down to 
fi = 3 GeV. The number of massless quarks for region B is set to four, and the heavy 
quark be should identified with the bottom quark, i.e., we have n¥°^ = 5. 
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^ (GeV) 


^/ 


as 




„MOM 












2 loop 


3 loop 




2 loop 


3 loop 


2 loop 


3 loop 


91.19 


5 


0.1180 


0.1180 


5 


0.1328 


0.1332 


0.1328 


0.1332 


50 


5 


0.1298 


0.1298 


5 


0.1480 


0.1486 


0.1480 


0.1486 


10 


5 


0.1779 


0.1781 


5 


0.2177 


0.2198 


0.2160 


0.2192 


4 


4 


0.2288 


0.2287 


5 


0.3074 


0.3148 


0.2988 


0.3096 


3 


4 


0.2536 


0.2538 


5 


0.3556 


0.3682 


0.3404 


0.3574 



Table 3: and a^^^ for various values of /i from region B. As input, as{Mz) = 0.118 
is used, which is transformed with the help of Eq. (jl]) to a^'^^{Mz)- The running to 
different values of /i is achieved with the help of the appropriate (3 function. a^^^(Q;s(/i)) 
is obtained from ^^(/i) using Eq. (jl]). 



(GeV) 




(nr 

a; 




„MOM 












2 loop 


3 loop 




2 loop 


3 loop 


2 loop 


3 loop 


91.19 


5 


0.1180 


0.1180 


5 


0.1180 


0.1181 


0.1180 


0.1181 


50 


5 


0.1298 


0.1298 


5 


0.1298 


0.1300 


0.1298 


0.1300 


10 


5 


0.1779 


0.1781 


5 


0.1797 


0.1804 


0.1798 


0.1807 


4 


4 


0.2288 


0.2287 


5 


0.2350 


0.2376 


0.2352 


0.2385 


3 


4 


0.2536 


0.2538 


5 


0.2612 


0.2652 


0.2615 


0.2667 



Table 4: as and a^^^ for various values of fi from region B. As input, as{Mz) = 0.118 
is used, which is transformed with the help of Eq. ( fT3l) to q;^°^(M^). The running to 
different values of is achieved with the help of the appropriate (3 function. a^^^{as{fi)) 
is obtained from ^^(/i) using Eq. (IT^ . 



In Tab. [T], we compare the values for for some selected fi values from region A in 
the MS and MOM schemes. For all numbers, we choose ai^\Mz) as the starting point, 
transform at /i = Mz to the MOM scheme and use the corresponding renormalization 
group equation to arrive at the desired n values. The same comparison for the case of the 
MOM scheme is shown Tab. |2l In both tables, we show in the last two columns the results 
of (Tab. m> and a™(/i) (Tab. [2]) as obtained from the MS-evolved value a,(/i) 

using Eqs. (jlj) and f|T3|) . respectively. The corresponding results for region B are shown 
in Tabs. [3]and|4] (where, of course, the values given in Eq. f|T9|) have been used). 

All four tables show good agreement between the values of the MOM coupling constant 
obtained via direct integration of the (quark- mass-dependent) MOM f3 function and with 
the help of the (simpler) conversion from the MS scheme. For Tabs. [1] and [3] the agreement 
is even getting better after taking into account the three-loop corrections. Note that the 
results in the case of the MOM scheme become slightly worse after switching on the three- 
loop terms, as can be seen in Tabs. [2] and HI The reason for this can be seen by comparing 
Eq. (|T6ll with Eq. (|TOl) . The former has (by construction) vanishing order corrections 
and a two-loop coefficient which is smaller by a factor of six. However, the three-loop 
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11.5 ^ 




o o o ooooooo 

o o o ooooooo 

|Li (GeV) ^ 

Figure 3: l/os as a function of /i. The (blue) upper solid line containing a step for 
/i = 2Mt corresponds to the MS result and the (red) lower solid line to the result in the 
MOM scheme. The (pink) dashed line represents l/a^ in the MOM scheme. The (black) 
dotted lines lying on top of the MOM and MOM result correspond to the results obtained 
from MS value of as using the conversion formulae (jlj) and ( JT3l) . respectively. For all 
results the three-loop expressions have been used and n^*-*^ = 6 has been chosen. 



term is only reduced by a factor of three and thus has bigger relative influence. Still, the 
difference between a^°^(/i) and ^^^^(^^(/i)) for = 3 GeV is about a factor of ten less 
than the current best value obtained, e.g., from hadronic r decay (see, e.g., Ref. [30]). 

In Figs. [3] and HJ the results of Tabs. [THU are shown in graphical from. In particular, 
we plot the inverse strong coupling as a function of both for the MS, MOM and MOM 
schemes, where in all cases the three-loop approximation is used for the running and the 
conversion between the schemes. We again choose ah\Mz) as the input quantity and 
convert at this scale to the other two schemes. The evolution of the MS coupling to lower 
/i values is shown by the (upper) solid lines with a step at the values for /if = 2Mf and 
lib = 2Mf,, respectively. Numerically very close is the dashed curve in the MOM scheme, 
which is expected from the above discussion. The lower solid line represents the result 
in the MOM scheme. Both for the MOM and MOM results, the conversion is performed 
for fi = Mz, and the running to other values of fi is achieved using the corresponding j3 
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10^ 



|Li (GeV) 

Figure 4: The same coding as in Fig. [3] has been adapted, except that ny^^ = 5 has been 
chosen. 



function. The dotted hues on top of the MOM and MOM curves represent the results 
where the transformation from the MS values is performed just at the considered value 
of n- 

5 Conclusions 

We have computed the three-loop corrections to the jS function of QCD with one heavy 
and 111 massless quarks in a momentum subtraction scheme (MOM). In our three-loop 
calculation, we do not consider the diagrams involving two different quark masses. Al- 
though there are only a few diagrams of this type, their evaluation is significantly more 
difficult. 

We have shown that our results describe the MOM coupling constant evolution in well- 
defined kinematical regions with three-loop accuracy. Moreover, the numerical analysis 
of our results has clearly demonstrated the full equivalence of the schemes with explicitly 
built-in decoupling (MOM and MOM) to the standard MS scheme, which, as is well- 
known, does not have such a property. 

From the more technical point of view, it has been shown that one can use the MS 
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scheme evolution along with simple conversion relations (derived for the regions either 
significantly above or below the heavy-quark threshold) to relate the values of the MOM 
scheme coupling constant from both regions. 

Finally, we believe that our analysis should help to remove the last traces of doubt 
about the usefulness of MS-like schemes, which formally do not obey the the Applelquist- 
Carazzone theorem [27j, for a unified description of mass effects in a broad region of /x^ 
values, from far below heavy-quark thresholds to well beyond them. 
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